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Based on the proper-time renormalization group approach, the scalar and the quark number 
susceptibiUties in the vicinity of possible critical end points of the hadronic phase diagram are 
investigated in the two-flavor quark-meson model. After discussing the quark-mass dependence of 
the location of such points, the critical behavior of the in-medium meson masses and quark number 
density are calculated. The universality classes of the end points are determined by calculating 
the critical exponents of the susceptibilities. In order to numerically estimate the influence of 
fluctuations we compare all quantities with results from a mean-field approximation. It is concluded 
that the region in the phase diagram where the susceptibilities are enhanced is more compressed 
around the critical end point if fluctuations are included. 
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I. INTRODUCTION 

Theoretical studies reveal an increasing richness in the 
structure of the phase diagram of strongly interacting 
matter. At high temperature T and/or quark chemi- 
cal potential fj, it is expected that the system undergoes 
a phase transition from the ordinary hadronic phase to 
a chirally restored and dcconfincd quark gluon plasma 
(QGP) 1113. Other phases such as a two-flavor (2SC) 
and color-flavor locked (CFL) color-superconducting 
phase are predicted at high densities and small tem- 
peratures. For small chemical potential and high tem- 
perature recent lattice calculations suggest the existence 
of mesonic bound states even above the deconfinement 
phase transition temperature 0, indicating a strongly 
coupled QGP. 

Based on model calculations 0,[EIE0 as well as lattice 
QCD simulations [111 113 , the existence of a critical end- 
point (CEP) )in the phase diagram is suggested. This is 
the endpoint of a first-order transition line in the (T, fi)- 
plane and is a genuine singularity of the QCD free energy. 
Here the phase transition is of second order, belonging to 
the three-dimensional Ising universality class. Its precise 
location, which is e.g. highly sensitive to the value of the 
strange-quark mass, is not known at present. It might be 
accessible with current and future experimental facilities 
and its observable implications in relativistic heavy-ion 
experiments such as event-by-event fluctuations of suit- 
able observables are intensively discussed ^ . 

From lattice studies it is known that the dynamics of 
the transition along the temperature axis is strongly af- 
fected by the presence of light quarks. The transition 
temperature is lowered substantially from its value in 
the pure gauge limit of infinitely heavy quarks, mostly 
due to the nearly massless up- and down quarks. This 
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indicates a prominent role of the (nearly exact) chiral 
SU{2) L X SU{2) fj symmetry, which is spontaneously bro- 
ken to SU{2)y at small temperatures and baryon densi- 
ties. Indeed, a chiral phase transition is observed numer- 
ically whose transition temperature coincides with that 
for the deconfinement transition. In the limit of van- 
ishing up- and down quark masses and infinite strange 
quark mass, the chiral phase transition is likely to be of 
second-order at vanishing /i and the static critical be- 
havior is expected to fall in the universality class of the 
Heisenberg 0(4) model in three dimensions [T^ . 

Lattice calculation of the QCD phase transition with 
finite ^ is much more difficult. Due to the Fermion sign 
problem, simulations of similar precision as at /i = 
are still lacking. From direct numerical evaluation of the 
QCDpartition function and its quantum-statistical anal- 
ysis [3, 13 or fi'om a Taylor expansion of the pressure 
around ^ = there is, however, evidence for a CEP at 
finite fi. At this point both the chiral- and the quark 
number susceptibility diverge. The existence of a CEP 
in the (T, ^)-plane implies that massless two-flavor QCD 
has a tricritical point (TCP) at which the 0(4) line of 
critical points ends. For larger /x's the transition is then 
of first-order. 

In an attempt to interprete the physical content of the 
lattice results there is a variety of model studies of the 
CEP and the critical region around it jl^, El, EE EE 
IT^ ITsl . In future searches for the CEP in heavy- ion 
reactions the size of the critical region is especially im- 
portant E3- Most of these studies rely in one way or 
another on a mean-field description of the phase tran- 
sition. As is well known, mean-field theory very often 
fails to give the correct description of phase transitions 
and fluctuations have to be included. An efficient way 
to describe critical phenomena beyond mean-field theory 
is the renormalization group (RG) method. It can be 
used to describe the universal and non-universal aspects 
of second-order as well as first-order phase transitions. In 
the context of the phase diagram of strongly interacting 
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matter the RG-method has been appUcd in |20j using a 
two-flavor quark-meson model, which captures essential 
chiral aspects of QCD. 

Extending our previous analysis, the present work pri- 
marily focuses is the question how large the critical re- 
gions around the TCP and CEP might be. The size of 
a critical region is defined as the one where mean-field 
theory breaks down and non-trivial critical exponents 
emerge. A well-known criterion is the Ginzburg criterion 
which estimates the size by some unknown coupling coef- 
ficients. It is based on an expansion of the singular part of 
the free energy for a second-order phase transition. Since 
the expansion coefficients, appearing in this criterion are 
not known for strong interactions the Ginzburg criterion 
is of limited use in the present context. Also universality 
arguments are not helpful, if the underlying microscopic 
dynamics is not well determined. This can be seen, for 
example, in the A transition of liquid helium He* and the 
superconducting transition of metals. Both transitions 
belong to the same universality class of the 0{2) spin 
model but their critical regions defined by their corre- 
sponding Ginzburg-Levanyuk temperature tgl, differ by 
several orders of magnitude. For hadronic matter it is 
expected that the critical region of the CEP is small ^3 ■ 

The paper is organized as follows. In the next Sec. Illll 
a mean-field analysis of the quark-meson model for two 
quark fiavors as an effective realization of the low-energy 
sector of QCD is presented. After the derivation of the 
grand canonical potential for the model, we calculate the 
phase diagram and localize the CEP. The behavior of the 
scalar- and pseudoscalar meson masses, the quark num- 
ber density and the scalar and quark number susceptibil- 
ities near the CEP are investigated. In order to obtain 
a deeper understanding of the shape of the critical re- 
gion around the CEP, critical exponents are calculated. 
In Sec. IIVI we repeat the calculation using the proper- 
time rcnormalization group approach (PTRG) in order 
to assess the infiucncc of fiuctuations. Finally, before 
summarizing in Sec. the size of the critical regions 
around the CEP and TCP are determined and compared 
with the mean-field results. 



component field 0. Without the explicit symmetry break- 
ing term c in the potential the Lagrangian is invariant 
under global chiral SU{2)l x SU{2)ii rotations. 

The four parameters of the model are chosen such that 
the chiral symmetry is spontaneously broken in the vac- 
uum and the a field develops a finite expectation value 
(cr) = /tt, where /jr = 93 MeV is the pion decay con- 
stant. Due to the pseudoscalar character of the pions the 
corresponding expectation values vanish (tt) = 0. 

The Yukawa coupling constant g is usually fixed by the 
constituent quark mass in the vacuum g = Mqj f^^. Us- 
ing the partially conserved axial vector current (PCAC) 
relation, the explicit symmetry breaking parameter c is 
determined by c = A/^/^, where M^^ is the pion mass. 
The quartic coupling constant A is given by the sigma 
mass M„ via the relation A = ^piM'^ — Ml). Finally, the 

parameter is found by minimizing the potential in ra- 
dial direction, yielding — o-2 _ c/(A(t). For the ground 
state where (cr) = J-^ this expression can be rewritten as 
= — M^/X . It is positive in the Nambu-Goldstone 
phase. 



III. MEAN-FIELD APPROXIMATION 

A. The thermodynamic potential 

In a spatially uniform system the grand canonical po- 
tential (per volume) is a function of the tempera- 
ture T and of the quark chemical potential fj,. We 
confine ourselves to the SU (2) /-symmetric case and set 
M ^ Mil = Md- This is a good approximation since fiavor 
mixing between the u- and d-quark in the vector channel 
is very small. 

The grand canonical potential is obtained as the loga- 
rithm of the partition function, which is given by a path- 
integral over the meson and quark/antiquark fields 



il/T 
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II. THE QUARK-MESON MODEL 

The Lagrangian of the linear quark-meson model for 
Nf = 2 light quarks q — (m, d) and Nc = 3 color degrees 
of freedom reads 

C = q{ip- g{(T + i'y5Tn))q 

+ ^{d^ad^^<j + d^TTd^Ti) - U{a, 7?) (1) 

where the purely mesonic potential is defined as 

C/((7,Tf) = ^(a2+7r2_^2)2_^^ _ (2) 

The isoscalar-scalar a field and the three isovector- 
pseudoscalar pion fields tt together form a chiral 4- 



To start with we evaluate the partition function in the 
mean- field approximation similar to Ref. |2lj . Thus we 
replace the meson fields in Eq. by their expectation 
values in the action neglecting the quantum and ther- 
mal fluctuations of the mesons. The quarks/antiquarks 
are retained as quantum fields. The integration over the 
fermions yields a determinant which can be calculated by 
standard methods (see e.g. Ref. H^l). This generates an 
effective potential for the mesons. Finally, one obtains 

n{T,y) = -Z^ = [/((a),(^)) + r!,-,(T,/i) (4) 
with the quark/antiquark contribution 



n^qiT, = VqT Jj^ {ln(l - n,(T, fi)) 
+ \n{l-ngiT,^l))} 



(5) 
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where 



• exp((i?, - fi)/T) 



(6) 



resp. nq{T, /i) are the usual quark/antiquark occupation 
numbers with nq(T,fi) = nq{T, —fi). The single-particle 

energy is given by Eq — ^Jk"^ + and the effective 
constituent quark mass by Mq = g{o'). The number of 
internal quark degrees of freedom is denoted by i>q — 
2NcNf = 12. 

The divergent vacuum contribution which results from 
the negative energy states of the Dirac sea has been ne- 
glected in Eq. (0) assuming that the renormalization is 
already done in the vacuum. This omission of the vac- 
uum contribution is in contrast to a similar analysis in 
the Nambu-Jona-Lasinio model where the vacuum con- 
tribution in the corresponding grand canonical potential 
cannot be neglected due to the finite momentum cutoff 
regularization [2]] |. 

Furthermore, we remark here that the same 
quark/antiquark contribution of the grand canonical po- 
tential can also be obtained by an analytical integra- 
tion of the corresponding RG flow equation in mean- 
field approximation |20l| . In this simple approxima- 
tion the mesonic RG flow is neglected and the residual 
quark/antiquark flow is given by 



tanh( 



2T ' 



tanh(^ 



2T - 

~i^E~ "'^9'^^''^) ~nq(T,n)] 
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which, after the scale integration towards the infrared, 
is identical to Eq. (jSj) if the vacuum contribution is ne- 
glected. 

For the case of a massless free quark/antiquark gas 
Eq. ||SJ) can be integrated analytically with the result 
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In general, any phase transition is characterized by an 
order parameter (j) which is identified here with the ex- 
pectation value of the sigma field. Its behavior is deter- 
mined by the corresponding classical equation of motion. 
It is obtained by minimizing the potential il in radial 
cr-direction. This results in the gap equation 

— = = X{a){{af -v')-c + g{qq) (9) 

where the scalar quark density (qq) has been introduced 
d^k 1 



{qq) = g{u)Vq 



(27r)3 



{n,(r,M)+ng(T,/.)} (10) 
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The quark density represents the source term in the 
equation of motion and can also be calculated analyti- 
cally for vanishing masses. In this limit, one obtains for 
the quark density 



M„^0 M„ 



lim 



12 



47r2 



(11) 



B. The phase diagram 

The solution of the gap equation © determines the be- 
havior of the order parameter as a function of the temper- 
ature and the quark chemical potential and thus allows 
to study the phase structure of the underlying two-flavor 
quark-meson model. 

In the vacuum we fix the model parameters to A/^r = 
138 MeV, Ma = 600 MeV, 
MeV which resuh in c 1.77 • 10^ MeV^, v 
A - 19.7 and g - 3.2. 
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FIG. 1; The phase diagram of the linear quark-meson model 
with two constituent quarks in mean-field approximation for 
three different fits of the pion mass. From top to bottom: 
= 138, 100, 69 MeV. 

For a physical pion mass (M^- = 138 MeV) the model 
exhibits a smooth crossover on the temperature axis 
and a first-order chiral phase transition on the density 
axis. For increasing temperatures this first-order transi- 
tion line persists up to a critical endpoint (CEP) where 
the chiral transition becomes second-order. Along the 
line of a first-order phase transition the thermodynamic 
potential has two minima of equal depth. These min- 
ima are separated by a finite potential barrier making 
the potential non-convex, which is typical for a mean- 
field approximation. The height of the barrier is largest 
at zero temperature and finite quark chemical potential 
and decreases towards higher temperature. At the criti- 
cal endpoint the barrier disappears (no latent heat) and 
the potential flattens. For temperatures above the CEP 
the transition is washed out and a smooth crossover takes 
place. With the parameters chosen above the location 
of the CEP is found to be at = 91.4 MeV, ^ic = 
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223 McV. In Fig. the phase diagram in the complete 
(T, /i)-plane is shown. 

To study the influence of explicit chiral symmetry 
breaking we have varied the model parameters keeping 
the Yukawa coupling g fixed. The resulting phase bound- 
aries are also shown in Fig. [T] for A/^r = 100 MeV and 
A/jr = 69 MeV (half of the physical pion mass). Re- 
ducing the pion mass while keeping the Yukawa coupling 
fixed, the CEP moves towards the T-axis. This is also 
the case if the Yukawa coupling is increased and the pion 
mass is kept fixed. Already for the pion mass M^r = 69 
MeV the CEP disappears and chiral symmetry is restored 
via a first-order transition for all temperatures and quark 
chemical potentials. As a consequence this model docs 
not have a tricritical point in the chiral limit (see also 
App.C in |23l|'). This is in contradiction to universality 
arguments as well as lattice QCD simulations. At vanish- 
ing quark chemical potential, the effective theory for the 
chiral order parameter is the same as for the 0(4) model, 
which has a second-order phase transition . It is fur- 
ther expected that the static critical behavior falls into 
the universality class of the 0(4)-symmetric Hciscnberg 
model in three dimensions [l^ . 

But, as already stated in Ref. [23|, within the mean- 
field approximation the order of the phase transition in 
the chiral limit of the linear quark-meson model strongly 
depends on the values for the model parameters. The way 
how to extrapolate towards the chiral limit is not unique. 
As wc shall see below, the mean-field approximation itself 
is also questionable [2^ . 



> 

0) 




20 40 60 80 100 120 140 160 180 
T [MeV] 

FIG. 2; Tfie pion and sigma masses as function of the tem- 
perature for three different quark chemical potentials (/i = 0, 
Ai = ~ 223 and ^ = 300 MeV). In the vicinity of the CEP 
(at /ic and Tc ~ 91 MeV) the sigma mass vanishes. 



C. In- medium meson masses 

To gain further physical insight into the critical behav- 
ior of the model we have studied the 'in-medium' meson 
masses M^{T, n) and M„{T, /i), which encode the critical 
fluctuations of the matter. In Fig.|21the temperature de- 
pendence of the meson masses for three different chemical 
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FIG. 3: Similar to Fig.|5| The pion and sigma masses as func- 
tion of the chemical potential for three different temperatures 
(T = 0, T = Tc ~ 91 and T = 150 MeV). 



potentials are shown. The sigma mass always decreases 
with temperature in the hadronic (chirally broken) phase 
and increases again at high temperature. The pion mass 
does not vary much below the transition but increases at 
high temperatures similar to the sigma mass indicating 
the restoration of chiral symmetry. The increase of the 
sigma and pion mass at high temperature is driven by 
the asymptotic behavior of the scalar quark density (see 
Eq. (|ll|l ') because chiral symmetry restoration enforces a 
decrease of the constituent quark mass. At large temper- 
ature the meson masses degenerate and increase linearly 
with T . A similar behavior is seen for the correspond- 
ing screening masses on the lattice at /i = |23|. For 
very high T ^ Tc the meson screening masses approach 
27rT because they are essentially controlled by the low- 
est Matsubara frequency of the quarks [22| . But for 
temperatures slightly above Tc the screening masses in 
the scalar and pseudoscalar channels are both still below 
27rT, which indicates a rather strong q-q interaction in 
the high-T phase. 

Fig. O shows the chemical potential dependence of the 
meson masses for three different temperatures analogous 
to Fig. 121 For small temperatures a discontinuity is visi- 
ble in the evolution of the masses, signaling a first-order 
phase transition. At T = and small /x the masses are 
constant and equal to the vacuum masses which can also 
be seen analytically. At the CEP the phase transition 
is of second order. Second-order phase transitions are 
characterized by long-wavelength fluctuations of the or- 
der parameter. Since the order parameter is proportional 
to the scalar cr-field, the corresponding scalar sigma mass 
must vanish at this point. In the vicinity of this point the 
sigma mass drops below the pion mass and the potential 
flattens in the radial direction. This behavior can clearly 
be seen in both Figs. [3 and 01 Since chiral symmetry is 
still explicitly broken the pseudoscalar mass stays always 
finite. For temperatures above the chiral transition the 
sigma mass drops below the pion mass and increases with 
the chemical potential (cf. Fig. |21and|3I). Close to /j,c for 
/i = 300 MeV in Fig.[2|the sigma mass drops rapidly with 
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temperature and jumps at the chiral transition. 



D. Susceptibilities 

In order to find a bound to the size of the critical region 
around the CEP wc calculate the quark number suscep- 
tibility Xq s-nd its critical behavior. In general, the quark 
number susceptibility is the response of the quark num- 
ber density n(T, ^) to an infinitesimal variation of the 
quark chemical potential 



(12) 



Before we study the thermodynamic properties of the 
susceptibility, we start with an analysis of the (total) 
quark number density n{T,ii). It is defined by a deriva- 
tive w.r.t. the quark chemical potential of the grand 
canonical potential (@J 



(27r)3 



{nq{T,^) ~ni,{T,^i)} . (13) 



For T ^ and /i ^ respectively, the analytical 
results for the quark number density are 

n(0, m) = - Mlf'^ 0{^, - M,), (14) 



n(T,0) = . 



(15) 



For a massless free quark gas one obtains for the density 

,2 " 



lim n(T,ri) ~ —n 



(16) 
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FIG. 4: The quark number density n in mean-field theory 
as a function of the chemical potential /i around the CEP for 
three difi'erent temperatures. The temperatures are Tc ~ 91 
MeV and r = Tc ± 5 MeV. 

In Fig. 21 the quark number density in mean-field ap- 
proximation is shown for three different temperatures 



around the CEP. For T = Tc the slope tends to infinity 
at /X = /Lie which will yield a diverging susceptibility. For 
temperatures below the critical one the system under- 
goes a first-order phase transition and the quark number 
density jumps. The slope of the quark density in the 
chirally symmetric phase is almost constant for all tem- 
peratures around the CEP. For temperatures above the 
CEP the discontinuity vanishes at the transition and the 
density changes gradually due to the smooth crossover. 
This produces a finite height of the quark number suscep- 
tibility Xq which is calculated via Eq. H12|l . The resulting 
susceptibilities are shown in Fig. [S] 
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FIG. 5: The quark number susceptibility in mean-field ap- 
proximation as a function of the chemical potential around 
the CEP for three temperatures. The temperatures are the 
same as in Fig.|l] (See text for details). 

For /.J = the susceptibility Xg(r, 0) increases as 
the quark mass Mq{T) decreases and reaches NjT"^ at 
Mq{T) = for three colors. Xq diverges only at the 
CEP and is finite everywhere else. The height of Xq de- 
creases for decreasing chemical potentials above the CEP 
towards the T-axis. For temperatures below the CEP Xq 
is discontinuous and jumps across the first-order transi- 
tion line. 

It is instructive to note that the quark number suscep- 
tibility is proportional to the isothermal compressibility 
kt via the relation 



i^T = Xq{T, n)/n'^{T, n) 



(17) 



In the vicinity of the critical point the quark number den- 
sity n is always finite but the susceptibility becomes large. 
This indicates that the system is easy to compress around 
the critical point and suggests that the interaction be- 
tween all constituents becomes very attractive |30f . 

It is also possible to express the quark number suscepti- 
bility as an integral over the quark/antiquark occupation 
numbers JH)) 

Xq{T,^i) = '^J-^-^{nq{T,^){l-nq{T,^i)) + 

nq{T,ti){l-nq{T,fi))} (18) 
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with the corresponding limits 



Xg(0,Ai) = ^MVM^^^^(/^-^^9) ' (19) 
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FIG. 6: The contour regions for two different ratios of 
the quark number susceptibiUties Rg = Xq{T , ^j) / 'x!q°° 0^ , 
{Rq — 3 and 5) in mean-field approximation around the CEP. 



As a function of temperature or quark chemical potential 
the maximum of Xa should coincide with the most rapid 
change in the chiral order parameter. One easily verifies 
that the scalar susceptibility is related to the sigma mass 
via ~ 

Similar equations can also be formulated for the pseu- 
doscalar susceptibility Xn, which is related to the order 
parameter via a Ward identity 

Xtt = {qq)/mq . (25) 

In the chiral limit the divergence of the pseudoscalar 
(transverse) susceptibility Xir in the broken phase signals 
the appearance of massless (Goldstone) modes. 

In Fig. [3 the contour region of the normalized scalar 
susceptibility 

is shown for four fixed ratios around the CEP. Again we 
observe an elongated critical region in the phase diagram, 
where Xa is enhanced in the direction parallel to the first- 
order transition line. The deeper reason for this feature 
can be understood by a study of the critical exponent of 
the susceptibility at the critical point which will be done 
in the next Section. 



For an estimate of the critical region around the CEP 
we have calculated the dimensionless ratio of the suscep- 
tibility and its massless free quark gas limit (cf. Eq. (|21|) ) 



XqiT,fl) 



(22) 



Fig. El shows a contour plot for two fixed ratios Rq in 
the phase diagram near the CEP. A constituent quark 
mass of roughly 300 McV has been chosen. The region 
of enhanced Xq is elongated in the direction parallel to 
the first-order transition line. 

In order to compare this behavior of the quark num- 
ber susceptibility around the CEP we repeat the previ- 
ous analysis with the scalar susceptibility. In general, 
static susceptibilities are obtained of the dynamic re- 
sponse function Xabi'^^O) in the static and long wave- 
length limit 



Xab 



1 d^n 

Vdadb 



lim Xab (0,g) 

9—0 



(23) 



where a, b denote external fields [23. The scalar suscep- 
tibility XiT corresponds to the zero-momentum projection 
of the scalar propagator, which encodes all fluctuations of 
the order parameter. The scalar susceptibility is related 
to the order parameter by 



£>{qq) 
drua 



(24) 
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FIG. 7: Same as Fig. |S|for four different ratios of the scalar 
susceptibilities R, = x<t (T, /x)/x<t (0, 0) {Rs = 10,15,20,25). 



E. Critical exponents 

Since the quark-meson model in mean-field approxi- 
mation does not exhibit a tricritical point in the chiral 
limit we focus on the critical behavior of the susceptibil- 
ity in the vicinity of the CEP. At that point the quark 
number susceptibility diverges with a certain critical ex- 
ponent. But a crucial observation is the following: In 
general, the form of this divergence depends on the route 
by which one approaches the critical point ■ For the 
path asymptotically parallel to the first-order transition 
line the divergence of the quark number susceptibility 
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scales with an exponent 7,. In mean- field approximation 
one expects 7^ = 1 for this path. For any other path, 
not parallel to the first-order line, the divergence scales 
with the exponent e = 1 — 1/S. Thus, in mean-field ap- 
proximation e = 2/3 because (5 = 3 and therefore 7g > e. 
This is the reason for the elongated shape of the critical 
region in the phase diagram (cf. Figs. |3 and 01 and why 
Xq is enhanced in the direction parallel to the first-order 
transition line. 




log (ln-n, |/[MeV]) 

FIG. 8: Tlie logarithm of the quark number susceptibility Xg 
as a function of In — /ic| at the CEP in mean-field approxi- 
mation. The solid line is a linear fit. 

In order to confirm this behavior we have calculated 
the critical exponent of the quark number susceptibility 
Xq numerically. We have used a path parallel to the fi- 
axis in the (T, /i)-plane from lower fi towards the critical 
fic ^ 223 MeV at fixed temperature ~ 91.4 MeV. 
Fig. |S1 shows the logarithm of Xq as a function of fi close 
to the CEP for a fixed constituent quark mass. Using a 
linear logarithmic fit we obtain 

InXg = - Mc| + r , (27) 

where the term r is independent of /i. We observe scal- 
ing over several orders of magnitude (only two orders 
are shown in the Figure) and obtain e = 0.64 ± 0.02, 
which is in good agreement with the mean-field predic- 
tion e = 2/3. The scaling close to the critical point is also 
demonstrated in Fig. |H| The onset of scaling is around 
In |yu — yUc| < —1 which is not shown in the Figure. 

Furthermore, other critical exponents can also be cal- 
culated by using susceptibilities. Considering the ratio 
R of the scalar and pseudoscalar susceptibility at finite 
temperature it is possible to calculate the critical expo- 
nent S [3^ . i? is a function of the reduced temperature t 
and the bare quark mass ruq 

R{t, mq) = x^'/Xa' = M^/M'^ ■ (28) 

At the critical point {t = 0) the order parameter scales 
as (qq) ~ mj^^. Using Eqs. TI^i and (^5)1 one obtains 

i?(0,m,) = l/,5, (29) 



which, at the critical point, is independent of rUq and 
yields (5 = 3 in mean-field approximation. In order to 
obtain this value it is necessary to determine the critical 
trajectory exactly at i = because there are two different 
limiting values for R due to spontaneous chiral symmetry 
breaking and its restoration. In the chiral limit the pions 
are Goldstone bosons in the broken phase (t > 0) yielding 
R{t > 0, 0) = and in the symmetric phase {t < 0) 
and Mo- become degenerate at large temperatures giving 
R{t < 0,0) = 1. 



IV. PROPER-TIME RG APPROACH 

In this section we study the phase diagram of the 
quark-meson model within the PTRG approach in the 
chiral limit as well as for realistic pion masses. In or- 
der to study the influence of fluctuations we will repeat 
the calculations of all previous quantities near the critical 
endpoint in the phase diagram. 

As we have seen, the mean-field approximation fails 
to properly describe the expected critical physics in the 
chiral limit (at least for the parameter set chosen). We 
will show, that this is remedied in the RG approach and 
a second-order transition in the chiral limit at finite tem- 
perature and zero chemical potential is found. This tran- 
sition lies in the 0(4) universality class, as expected. For 
nonzero chemical potential the second-order transition 
ends in a TCP and becomes a smooth crossover for fi- 
nite quark (pion) masses with a CEP. Thus, by variation 
of the pion mass the relationship and the correlations 
between the TCP and the CEP's can be studied. In ad- 
dition, the influence of fluctuations on the susceptibilities 
and the critical region around the CEP can be assessed. 

In order to include fluctuations we use the proper-time 
renormalization group (PTRG) method [H Im HI HI 
113, 111] • In the vacuum, the PTRG flow equation for the 
scale-dependent effective action Tkl^] is governed by the 
selfconsistent equation 



(2) 

where F^. represents the full inverse propagator and is 
given by a second functional derivative w.r.t. the field 
components $. The smearing function is labeled by 
/fe(rfc^) and tr denotes a four-dimensional momentum 
integration and a trace over all given inner spaces (e.g. 
Dirac, color and/or flavor-space). Details and the gen- 
eralization to finite temperature and chemical potential 
can be found e.g. in Refs. [20ll3^ . 

For the quark-meson model 1^ the resulting flow equa- 
tion for the scale-dependent grand canonical potential 
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^lk{T,fi) is given by 



3 .1.1^^- 



-— cothf — ^ 



2NcNf 



with 

E„ = 



En 



the 



,(30) 



pion- 



E^ 



yWT2n^, the cr-meson 
and quark energies Eg ~ 



i/fc^ + 3^0^. The primed potential denotes the 0^- 
derivative of the potential, i.e., n'j, := d^lk/dcj)^ and cor- 
respondingly the higher derivatives. 

We solve the flow Eq. (|30|l by discretizing the generally 
unknown potential f2fe on a 0^-grid. For details concern- 
ing the numerical implementations we refer the interested 
reader to Refs. [20ll3^ and references therein. 

For the results presented below we have chosen an UV 
cutoff A = 500 MeV and have fixed the quartic coupling 
to A = 11.5 in order to reproduce a pion decay constant 
/tt ~ 87 MeV in the chiral limit. The running of the 
Yukawa coupling g is neglected here and fixed to g = 4.2 
in order to reproduce a constituent quark mass of the 
order of 370 MeV. Note, that this value of the quark 
mass is roughly 100 MeV larger than the one used in 
the previous work |20l |. As a consequence, the tricritical 
point in the phase diagram moves to larger temperatures 
and smaller chemical potentials. 



A. The phase diagram 

In the present work we have also generalized previ- 
ous results in the chiral limit |23 to finite pion masses. 
The resulting phase diagrams for the chiral limit and for 
physical pion masses ~ 130 MeV are both shown in 
Fig. In the chiral limit a second-order phase transition 
(dashed line) belonging to the 0(4)-universality class is 
found. For zero chemical potential we find a critical tem- 
perature Tc ^ 170 MeV for two massless quark flavors, 
in good agreement with lattice simulations |40l |. For in- 
creasing chemical potential the second-order transition 
line ends in a TCP (bullet). 

For our choice of parameters the location of the TCP 
is at r^* ~ 80 MeV and ^* = 268 MeV. For temperatures 
below the TCP the phase transition changes initially to a 
first-order transition (solid line). For temperatures below 
Ts ^ 10 MeV two phase transitions with a further tricrit- 
ical point (labeled as '2nd TCP' in the Figure) emerge 
as already described in Ref. [23|. A larger constituent 
quark mass pushes the location of the TCP towards the 
temperature axis and the location of the splitting point 
down towards the chemical potential axis. All qualitative 
properties of the two phase transitions below the splitting 
point survive. Only the area bounded by the transition 
lines is reduced for increasing quark masses (cf. phase 
diagram in Ref. [20l|). 

As the pions become massive the TCP turns into a 
CEP. For each value of the pion mass there is a cor- 
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FIG. 9: Two phase diagrams of the linear Nf = 2 quark- 
meson model: One for physical pion masses (right solid line 
which ends in the CEP) and another one for the chiral limit. 
Sohd lines denote first-order and dashed lines second-order 
transition lines (see text for details). 



responding CEP. In an extended, three-dimensional (T, 
/i, m,r) phase diagram these points arrange into a criti- 
cal line (the 'wing critical line'). It is expected that the 
static critical behavior of this line falls into the univer- 
sality class of the Ising model in three dimensions cor- 
responding to the one-component i^^-theory in 3D. For 
the coupling constants listed above Af^ ^ 123 MeV and 
Mq ^ 390 MeV in the vacuum and the location of the 
CEP is at - 61.5 MeV and ~ 313.1 MeV. A better 
fit to the physical pion mass can be accomplished by a 
fine tuning of the initial coupling constants. Compared 
to recent lattice and other model studies (cf. Fig. 6 in 
Ref. 0) the location of the TCP and consequently of the 
CEP is at lower temperatures due to omission of other 
degrees of freedom in the quark-meson model. 

For any finite value of the symmetry breaking param- 
eter c in the potential Q the critical 0(4)-line is im- 
mediately washed out and turns into a smooth crossover 
line which is not visible in the phase diagram. For tem- 
peratures below the CEP a first-order curved transition 
line is found which persists down to the ^-axis. In the 
chiral limit below the splitting point, the right second- 
order transition turns into a crossover for finite quark 
masses and is again not visible in the phase diagram any- 
more. Analogously, the second tricritical point, if it ex- 
ists, should turn into a critical point. Some remnants of 
this critical point can indeed be seen in the vacuum ex- 
pectation value and meson masses. But a detailed anal- 
ysis of this point is postponed to a future work. In the 
following we focus of the region around the TCP and 
CEP. 



B. In- medium meson masses 

The scale-dependent in-medium sigma and pion masses 
are determined by the curvature of the potential il at the 
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global minimum 0o via the relations 



<j>=(t>0 

(31) 
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FIG. 10: The pion and sigma masses as a function of tem- 
perature for two different quark chemical potentials = 
and = ^ij; ~ 313 MeV). Near the CEP (located at ^il and 
~ 61.5 MeV) M„ vanishes. 



In Fig. ^1 the meson masses are shown as a function 
of temperature for two different chemical potentials sim- 
ilar to Fig. 13 For /i = the sigma mass drops rela- 
tively quickly with temperature while the pion mass are 
almost constant in the chirally broken phase. This be- 
havior is very similar to a PTRG analysis with another 
smearing function and a truncation of the meson poten- 
tial (cf. Fig. 14 in Q). For = ^J: - 313 MeV and 
at r = the sigma mass is almost degenerate with the 
pion mass ~ 150 MeV (see Fig. llOf) . For increas- 
ing temperatures, the sigma mass jumps to almost its 
vacuum value as it crosses the first-order phase transi- 
tion. In the vicinity of the CEP around T ~ 60 MeV 
the sigma mass drops below the pion mass and vanishes 
at the critical point. At this point the potential in radial 
direction around the global minimum becomes flat. As in 
the mean-field approximation the slope of the sigma mass 
as function of temperature (or of the chemical potential) 
at the CEP tends to infinity. 

In Fig. ^2 the meson masses parallel to the /i-axis for 
three different temperature slices (T = 1, T = and 
T = 150 MeV) in the phase diagram are shown. For T = 
1 MeV the sigma mass is constant and slightly before the 
first-order transition a melting of the mass takes place. 

Due to the curved first-order transition line the criti- 
cal chemical potential at the CEP around T ^ &1 MeV 
coincides with the critical chemical potential at T = 1 
MeV. 

One important observation is that the sigma mass 
drops to zero around the CEP much faster if fluctuations 
are included. 
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FIG. 11: Similar to Fig. 1101 The pion and sigma masses as a 
function of chemical potential for three different temperatures 
(T = 1, T = T<f ~ 61 and T = 150 MeV). 



C. Susceptibilities 

In the vicinity of critical phenomena, fluctuations 
which arc neglected in a mean-field approximation be- 
come more and more important. In order to study the 
modifications, induced by the fluctuations, we have re- 
calculated the quark number density and susceptibilities 
within the RG method. It should be recapitulated that, 
in contrast to the mean-field approximation, this method 
yields a TCP in the phase diagram for the chiral limit. 

In the following we begin our investigations with the 
chiral limit. In Fig. ll2l thc density n is shown as function 
of /i for three different temperatures around the TCP. 



E 
c 
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260 262 264 266 268 270 272 274 
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FIG. 12: The quark number density around the TCP. The 
three temperatures are = 80.2 MeV and T = ± 5 MeV. 

For temperatures below the tricritical T*^ n jumps due 
to the first-order transition. Above the tricritical tem- 
perature there is a kink in the density due to the second- 
order nature of the transition. In the symmetric phase n 
grows almost linearly. At the TCP the slope of the den- 
sity when approached from the broken phase, diverges 
thus yielding a divergent susceptibility. This is demon- 
strated in Fig. El where the chemical potential depen- 
dence of the quark number susceptibility is shown for 
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three different fixed temperatures around the TCP. 



tures around the CEP, similar to Fig. is displayed. 




260 262 264 266 268 270 272 274 
H [MeV] 

FIG. 13: The quark number susceptibility around the TCP. 
The three temperatures are the same as in Fig. 1121 

In the chiral limit, the susceptibility always jumps 
across the first- or the second-order transition. In Fie . 1131 
the second-order transition jump of Xq is drawn with 
solid lines. Far below the chiral phase transition Xq is 
suppressed. In the restored phase Xq tends towards the 
value of the massless free quark gas, x,''"'', Eq. (|^ . 

When we now leave the chiral limit and investigate the 
phase diagram around the CEP, the critical behavior of 
the transition changes. In contrast to the behavior of the 
quark number density around the TCP it does not have 
any kink at finite quark masses (cf. Fig. I14|l. 
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FIG. 14: The quark number density around the CEP. The 
three temperatures are = Ql.^ MeV and T = T<f ± 5 MeV. 

For temperatures below the critical T^, the density 
jumps again at the phase transition due to its first-order 
character. But above the CEP the discontinuity vanishes 
at the transition and the system runs through a smooth 
crossover. For temperatures above T^, the slope of the 
density is always finite, except exactly at where it 
diverges yielding a divergent quark number susceptibil- 
ity at this point. Thus, the divergence of Xq survives 
even at finite quark masses. The finite slope results in 
a finite peak of Xq which is shown in Fig. 1151 where the 
quark number susceptibility for three different tempera- 




306 308 310 312 314 316 318 320 
H [MeV] 

FIG. 15: The quark number susceptibility around the CEP as 
a function of chemical potential for three fixed temperatures. 
The three temperatures are T^'^ = 61.5 MeV and T = T^'^ ± 5 
MeV. 



Below the chiral transition, the susceptibility is again 
suppressed and in the restored phase it tends towards 
X^'®® of a massless free quark gas. Only at the first-order 
transition there is a discontinuity. In contrast to the chi- 
ral limit, the quark susceptibility is a smooth continuous 
function for temperatures above the CEP. 

On the lattice a similar behavior of the quark suscep- 
tibility is seen At ^ = 0, Xq increases smoothly 
near the critical temperature. At finite chemical poten- 
tial it develops a pronounced dip but with increased error 
bars around the transition temperature. On the other 
side, a significant peak in the isovector channel is not 
seen on the lattice. Thus, this indicates that only one 
scalar degree of freedom becomes massless at this point. 
Also universality arguments predict that the susceptibil- 
ity diverges at both the TCP and CEP with certain but 
different critical exponents. At the TCP Xq should show 
a power-law behavior with a critical exponent 7g. Since 
the TCP belongs to a Gaussian fixed point mean-field 
exponents are expected. 

In the chiral limit the susceptibility Xq always has a dis- 
continuous jump across the critical 0(4)-transition line. 
Xq is larger below the 0(4)-line (in the chirally broken 
phase) than in the restored phase. At any finite /i the sus- 
ceptibility has a cusp at the critical 0(4)-transition line. 
The peak of this cusp on the critical line as approached 
from the broken phase becomes higher and higher as we 
increase /i. Finally, it diverges exactly at the TCP, while 
the potential and thus the pressure will stay finite. Only 
at /I = 0, the quark number susceptibility is a continuous 
function across the phase boundary. 

On the other side, as already mentioned, in the chi- 
rally symmetric phase, Xq tends towards the value of the 
massless free quark- gas susceptibility x^^'^'' . The devia- 
tions are caused by the residual meson interaction in this 
phase. 
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D. Critical exponents 

For finite quark masses it is expected that the univer- 
sality class of the critical points changes. In order to 
determine the universality class of the TCP and CEP 
we have also calculated the corresponding critical expo- 
nents of the quark number susceptibility. In determing 
the critical exponent we proceed in exactly the same way 
as described in the mean-field section. The calculation 
of the first and second derivative of the potential for the 
quark number density and susceptibility is again deter- 
mined numerically with an adaptive step size algorithm 
in order to minimize rounding and truncation errors. 

We again begin the analysis with the chiral limit. We 
have used a path parallel to the /i-axis in the (T, fi)- 
plane from lower ^ towards the TCP. Note, that in the 
(T, /x)-plane the first-order transition line and the critical 
0(4)-line at the TCP are asymptotically parallel and are 
not parallel to either the T- or /i-axis. 



2 
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FIG. 16: The logaritfim of the quark number susceptibility 
around the TCP as a function of the logarithm of the chemical 
potential in the vicinity of the TCP. The solid line is a linear 
fit. 

Using a linear logarithmic fit of Xq in the vicinity of 
the TCP we obtain the critical exponent = 0.53 which 
is shown in Fig. 1161 The region where the scaling starts 
is rather smaU |/i - | < 10~°-^ ~ 0.3 MeV. This expo- 
nent is consistent with mean-field theory which predicts 
a critical exponent 7, = 1/2. This is also expected due 
to the Gaussian fixed point structure of the TCP. 

Leaving the chiral limit, the universality class of the 
CEP is that of the three-dimensional Ising model. Using 
again a path parallel to the /x-axis from lower ^ towards 
the CEP we repeat the calculation of the critical expo- 
nent. In Fig. 1171 is shown versus the distance from the 
CEP. It is interesting to see, that there seems to be two 
different scaling regimes, which we indicated with two 
linear fit functions in the figure. 

The slope for the data points changes between 10~°-^ 
MeV < l/i - Mcl < 10°-^ MeV. We have fitted the data 
for III - /i^l < 10-°-5 MeV and > 10°-^ MeV separately 
and obtained the critical exponent e ~ 0.74 (solid line) 
for \fi - fi^J < 10-°-5 ~ 0.3 MeV. For |/i - fi^J > 10°-5 
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FIG. 17: Same as Fig.EHlfor the CEP. 



McV we also see a linear behavior for several orders of 
magnitude. If there is a proper scaling behavior of the 
susceptibility in this region, the slope would be consistent 
with an exponent of e ~ 0.77 (dashed line). This change 
of the exponents could be interpreted as a crossover of 
different universality classes ^3 • One possibility for this 
crossover phenomenon is that the CEP is affected by the 
TCP. In Ref. ^3 a similar crossover phenomenon be- 
tween different universality classes in the framework of 
the CJT potential for QCD in improved-ladder approx- 
imation is seen: As they approach the CEP for realis- 
tic quark masses (m^ 5 MeV), the critical exponent 
change gradually from those of the tricritical point to 
those of the 3D Ising model via those of the CEP in the 
mean-field approximation. In mean-field theory one ex- 
pects e = 7//3<5 = 2/3. They find a crossover from the 
nontrivial critical exponent e = 0.57 ± 0.01 to the mean- 
field exponent e = 0.68 ± 0.02. Contrary to this work, 
we find a crossover from e = 0.77 to e = 0.74. Thus, 
at the CEP the susceptibility diverges with the critical 
exponent e ~ 0.74. This exponent is consistent with the 
one of the 3D Ising universality class e = 0.78 and is def- 
initely different from the mean-field value e = 2/3. Note 
also, that the mean-field exponent of a bicritical point 
(CEP) are in general different from those of a tricritical 
point. In order to complete the analysis here we sum- 
marize the critical exponents of different approaches in 
Tab.Hl 





a 


/3 


7 


S 


V 


■n 


Ref. 


mean-field 





1/2 


1 


3 


1/2 





[43] 


e-exp. (5 loops) 


0.11 


0.327 


1.24 


4.80 


0.631 


0.04 


[44J 


rmm. sim. 


0.12 


0.31 


1.25 


5.20 


0.64 


0.06 


m 



TABLE I: Critical exponents for the 3D Ising model. 



It is also interesting that the value |/i — //J;| ~ 0.5 MeV 
where the slope of the data points changes to a mean-field 
scaling in Ref. [T^ is comparable to our value. 
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E. The critical region 

We are now able to compare the impact of the quan- 
tum and thermal fluctuations on the shape of the region 
around the critical points in the phase diagram. We re- 
peat the calculation of the contour plot of Figs. El and [7| 
in the framework of the RG approach. The result for the 
scalar susceptibility around the CEP for the same four 
ratios Rs as in the mean-field approximation is shown in 
Fig. Hi 




-0.05 0.05 0.1 



FIG. 18: The contour regions for four different ratios of the 
scalar susceptibilities X<t (T, m)/x<t (0, 0) [Rs = 10,15,20,25) 
around the CEP in the phase diagram. 

As in the mean-field case the region is elongated in 
the direction of the first-order transition line, but it is 
now much more compressed. For example, choosing a 
ratio of i?s = 15, the corresponding susceptibility covers 
an interval from —0.04 to 0.02 in the reduced chemical 
potential direction and an interval from —0.1 to 0.4 in 
the temperature direction. In the mean-field case the 
same ratio covers an interval from —0.8 to 0.1 in the fi 
direction and from —0.15 to 0.6 in the T direction. While 
the interval in the temperature direction is comparable in 
both cases, the effect in the chemical potential direction is 
dramatic. In the RG calculation the interval is shrunken 
by almost one order of magnitude, despite the fact that 
the corresponding critical exponents are quite similar! 

An similar result is obtained for the critical region of 
the quark number fluctuations. In order to compare the 
critical region around the CEP with the one around the 
TCP we show the quark number susceptibility in Fig. 1191 
in a larger sector of the phase diagram containing both 
critical points. 

The CEP is far away from the TCP at larger values 
of the chemical potential and smaller temperature as ex- 
pected. Due to the sharp transition lines in the chiral 
limit the critical region around the TCP is chop off in 
the chirally symmetric phase. Whether there is a robust 
effect of the TCP on the CEP as stated in Ref. [l3| is 
not seen in this work. But this important issue includ- 
ing the quark mass dependence of the susceptibilities is 
postponed to a further analysis. 
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FIG. 19: The contour regions for two different ratios of 
the quark number susceptibilities Rq = Xg{T, fJ-)/xt^°iT, /J^) 
around the CEP and TCP in the phase diagram. 

V. SUMMARY AND CONCLUSION 

Using a Wilsonian renormalization group approach, we 
have analyzed the phase diagram of hadronic matter in 
the two-flavor quark-meson model. This model captures 
essential features of QCD, such as the spontaneous break- 
ing of chiral symmetry in the vacuum and can therefore 
yield valuable insight into the critical behavior, associ- 
ated with chiral symmetry. Of special importance is 
the emergence of a CEP which is intensely discussed at 
present, in connection with fluctuation signals in heavy- 
ion collisions. Here the size of the critical region around 
the CEP is of special importance. Most studies of this 
issue, available in the literature, have been performed 
in the mean-field approximation, which neglects thermal 
and quantum fluctuations. These can be assessed, how- 
ever, in the RG approach which is able to correctly pre- 
dict critical exponents in the vicinity of critical points of 
the phase diagram. 

The main results of this work can be summarized as 
follows: 

1. From universality arguments it is expected that the 
quark-meson model (and most likely QCD) has a 
TCP in the chiral limit. For the parameter set cho- 
sen in this paper, a mean-field calculation is not 
able to find such a point, while the RG predicts 
its existence. The expected critical behavior is also 
reproduced. Due to the Gaussian fixed-point struc- 
ture at the TCP mean-field exponents are expected 
what we could verify. 

2. When effects of finite current quark masses (finite 
pion masses) are included, a CEP emerges in both 
the mean-field and the RG calculation. By analyz- 
ing the scalar- and quark number susceptibilities 
we find that the RG calculation yields non-trivial 
critical exponents, consistent with the expected 3D 
Ising universality class. Close to the CEP our ex- 
ponent is consistent with the Ising class but we also 
see a novel crossover phenomenon. 
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3. As a consequence of fluctuations, the size of the 
critical region around the CEP is substantially re- 
duced, as compared to mean-field results. This is 
particularly true in the //-direction, where a shrink- 
age by almost one order of magnitude is observed. 
This may have consequences for the experimental 
localization of the CEP in the phase diagram since 
it further complicates its detection through event- 
by-event fluctuations. 

The success of the RG approach in predicting the ex- 
pected critical behavior of the thermodynamics in the 
quark-meson model encourages us to pursue the issue of 
the existence of a CEP in the phase diagram of strongly 
interacting matter and its location in the (T, /t)-plane 
further. It is known from lattice studies that the strange 
quark plays a major role. We therefore intend to extend 
the analysis to three quark flavors . One of the draw- 



backs of the quark-meson model is the lack of explicit 
gluonic degrees of freedom, which are known to play a big 
role in the thermodynamics of QCD and are associated 
with conflnement aspects of the theory. One possibility 
to incorporate such effects is the coupling of the quark- 
meson model to Polyakov-loop model [ItI li^ , as pu t 
forward in . Work in this direction is in progress [Ho| • 
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